Interstellar Seeing. I. Superresolution Techniques Using Radio Scintillations 



J. M. Cordes 

Astronomy Department, 520 Space Sciences Building, Cornell University, Ithaca, NY 14853 

cordes@spacenet.tn.cornell.edu 

Submitted to Astrophysical Journal February 5, 2008 
ABSTRACT 

Interstellar scintillation can be used to probe transverse sizes of radio sources on 
scales inaccessible to the nominal resolution of any terrestrial telescope, e.g. ^ I0~^ 
arc sec. Methodology is presented that exploits this superresolution phenomenon for 
both single aperture and interferometer observations. The treatment applies to the 
saturated (strong-scattering) regime and holds for both thin screens and extended 
media. A general signal model for radio sources is presented, scintillated amplitude 
modulated noise, which applies to compact, incoherent synchrotron sources such as 
AGNs and gamma-ray burst sources and also to known, coherent sources such as 
masers and pulsars. The exact probability density function for measured intensities 
and interferometric visibilities is obtained by solving a general Fredholm problem. 
An approximate density function is also obtained by using the equivalent number of 
degrees of freedom in scintillation modulations. The scintillation modulation variance 
is presented, which includes the effects of source structure and time-bandwidth 
averaging in the signal processing. Two Bayesian methods are outlined for inferring 
the sizes of emission regions that use first order statistics of the intensity and 
visibility. Intensity cross-correlation methods for inferring source sizes are also given. 
Intensity interferometry in the the radio context is compared to the optical intensity 
interferometry of Hanbury-Brown and Twiss. 



1. Introduction 



Diffractive interstellar scintillation (DISS) is caused by multipath scattering of radio waves 
from small-scale density irregularities in the ionized interstellar medium. It is sensitive to intrinsic 
sizes of radiation sources in much the same way that optical scintillation from atmospheric 
turbulence is quenched for planets while strong for stars. However, interstellar scintillation differs 
from the atmospheric case in that it can resolve sources at angular resolutions much smaller 
than those achievable with available apertures, including the longest baselines used in very long 
baseline interferometry (VLBI), those using space antennas. Optical techniques such as intensity 
interferometry, speckle interferometry and adaptive optics typically only restore the telescope 
resolution to what it would be in the absence of any atmospheric turbulence. 
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We define the superresolution regime where the source is unresolved by terrestrial 
interferometers but is sufficiently extended to modify the DISS. Let Os,Oij and 9iso be the source 
size, interferometer fringe spacing, and isoplanatic DISS patch, respectively. By definition, 
two point sources separated by much less than the isoplanatic angle will show identical DISS. 
The isoplanatic angle Oiso ^/dOd, where A is the wavelength, Od is the size of the scattering 
("seeing") disk and d is the source-Earth distance. Using typical numbers {d = 1 kpc, 9d = I 
mas at an observing frequency of 1 GHz), 9iso ~ 0.4/Ltarc sec. For pulsars, whose light-cylinder 
radii, tlc = cP/2ir (P = spin period) arc smaller than 1 ^arc sec at typical distances, we have 
Giso ^ Os <^ Oij and 9iso ^ ^d- In this case, speckle methods can achieve far better resolution than 
the interferometer. By speckle methods, we mean observations that analyze differences in the DISS 
between source components, which are sensitive to the spatial separations of those components. 
Cornwell & Narayan (1993) have discussed particular superresolution techniques in the radio 
context. In optical astronomy, superresolution is not achievable because 9d ~ ^iso- However, 
the superresolution regime has been identified in optical laboratory applications (Charnotskii, 
Myakinin Sz Zavorotnyy 1990). 

The purpose of this paper is to provide a rigorous and general formulation of DISS that can be 
used in superresolution applications. Previous work has relied on DISS theories that are restricted 
to particular spatial geometries (e.g. thin screens) or to particular wavenumber spectra (e.g. a 
power law with the "Kolmogorov" slope.). In subsequent papers, we will apply the methods of 
this paper to scintillation observations and derive constraints on pulsar source sizes. 

Our treatment builds upon work published in both the optical and radio propagation 
literature and is unique in the following ways: (1) it interfaces empirical astrophysical constraints 
on source radiation fields with the scattering geometry and scattering strength appropriate for 
radio observations; (2) it calculates fluctuation statistics of realistically measureable quantities in 
single-aperture and interferometric observations, taking into account all sources of fluctuation; 
(3) we present exact calculations for the intensity probability density function (PDF) that take 
into account arbitrary source brightness distributions and arbitrary amounts of time-bandwidth 
averaging; and (4) our results can be applied to a wide range of media with arbitrary spatial 
extent and wavenumber spectrum. Gwinn et al. (1998) discuss issues that are very similar to 
those contained in this paper. Our treatment is more general than Gwinn et a/.'s because it is 
not limited to scattering media contained in thin screens. Also, our treatment of the probability 
distribution is not restricted to small source sizes (compared to the isoplanatic scale). Finally, 
our treatment includes the effects of intrinsic source fluctuations, which are assumed negligible by 
Gwinn et al. 

We restrict the analysis to strong (saturated) DISS. Another component — refractive 
interstellar scintillations (RISS) — also modulates source intensitities, but on time scales much 
longer than for DISS. RISS can be seen from larger angular diameter sources, by a factor of 1000, 
than can DISS. We are interested in modeling data spans much shorter than the characteristic 
RISS time scale, so we ignore RISS in our discussion. The cases we consider are what would be 
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called the "single speckle" regime in the optical literature. This corresponds to the case where 
the aperture size (either a single dish diameter or an interferometer baseline) is smaller than the 
diffraction length scale (e.g. the Fried scale). However, our formalism can be extended easily to 
include aperture averaging and multiple speckle cases. 

In §^ we briefly summarize previous work on use of DISS to resolve radio sources. In §^ 
we present a general signal model and in §^ give expressions for the modulation index of the 
interferometric visibility and single aperture intensity. We take into account time-bandwidth 
averaging and source extent. We calculate an approximate PDF of the visibility and intensity 
using the number of degrees of freedom in the scintillations. We also outline the exact calculation 
of the PDF. §^ presents a Bayesian inference method for the source size. The paper is summarized 
in §p. The details of our definitions and calculations are given in three Appendices. In Appendix 
we derive the statistics and scintillations of the scintillating amplitude modulated noise model. 
In Appendix ^ we derive second moments for the intensity and visiblity. In Appendix ^ we 
derive the PDFs for the scintillations, total intensity, and visiblity. 

Notation required in discussions of wave propagation through random media is necessarily 
copious. Apart from standard definitions for wavelength, frequency, speed of light, and wavenumber 
(A, z^, c and k = 27r/A, we list in Table 1 those symbols that are used throughout the paper. 

2. Previous Applications of DISS Superresolution 

DISS superresolution techniques have been used in several ways to probe source structure. 
DISS has been sought from various kinds of AGNs (Condon & Backer 1975; Armstrong, Spangler 
& Hardee 1977; Condon & Dennison 1978; Dennison &: Condon 1981), leading to bounds on the 
brightness temperature and limits on the Lorentz factors of bulk relativistic flow. Lovelace (1970) 
first suggested that DISS might resolve pulsar magnetospheres. Backer (1975) placed coarse limits 
on the sizes of pulsar magnetospheres by considering the fractional modulation of DISS for a few 
pulsars. Cordes, Weisberg &; Boriakoff (1983; hereafter CWB83) placed upper bounds on the 
source extent, and emission altitude, for two pulsars, finding upper bounds on emission altitudes 
of 0.5rLc and O.lrLc- Wolszczan & Cordes (1987) exploited a remarkable episode of multiple 
imaging of a pulsar by a refracting interstellar structure; they found that (a) different pulse 
components from a long-period pulsar showed nonidentical DISS, signifying source resolution; and 
(b) the implied emission altitude is comparable to tlci in contrast to other estimates, based on 
pulse widths and polarization, which suggest emission altitudes of only 1-10% of tlc (Blaskiewicz 
et al. 1991). Kuzmin (1992) and Smirnova, Shishov & Malofeev (1996) reported similar results 
on four additional long-period pulsars, again with implied emission radii ~ tlc- 

Recently, Gwinn et al. (1997,2000) have analyzed VLBI observations of the Vela pulsar using 
time and frequency resolutions that exploit the same spatial resolving power of DISS as do the 
single-dish observations used by others. Through estimation of the probability density function 
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(PDF) of the visibility function magnitude, they infer that DISS shows less modulation than 

expected from a point source and therefore conclude that the source must be extended. They 
estimate a transverse size ~ 500 km for the region responsible for the pulsed flux in a narrow 
range of pulse phase. 

In a second paper, we reassess the conclusions of Gwinn et al. (1997) by considering how 
time-frequency averaging in the signal processing affects inferences on source size from visibility 
fluctuations. Gwinn et al. (1999) have also considered time-frequency averaging and the effects 
of source fluctuations. The approaches differ and yield different conclusions about the importance 
of averaging and source noise as well as differing on estimates of the source size. 

Also recently, radio observations of gamma-ray burst (GRB) afterglows suggest that DISS 
occurs in the early stages and then is quenched as expanding synchrotron sources are first 
smaller than, and then exceed, the isoplanatic scale (Goodman 1997; Frail et al. 1997). If GRB 

synchrotron sources are incoherent sources, the angular size requirements for DISS to appear are 
severe. We defer to another paper a discussion of GRB scintillation. 



To account for all contributions to measureable quantities, we need a comprehensive statistical 
model for the received signal. The model presented here includes a source that is temporally 
incoherent but has arbitrary spatial coherence; diffractive interstellar scintillation from an 
arbitrary distribution of scattering material along the line of sight; and additive radiometer noise. 

We define e{r,t) to be the complex, narrowband, baseband scalar electric field that is 
explicitly or implicitly manipulated in radio astronomy systems. (For definitions, see Appendix A.) 
It is determined by the source emission mechanism and by propagation effects through intervening 
media, as well as by receiver and background sky noise. Wc consider source emission that has 
underlying Gaussian statistics and propagation effects from the turbulent, ionized interstellar 
medium (ISM). The field measured at position r is the superposition of scintillating source 
components and radiometer noise, n{r,t), 



All vectors are two-dimensional and transverse to the line of sight. The corresponding intensity is 



3. Signal Model for Scintillating Sources 




(1) 



I{r,t) = \e{r,t)\\ 



(2) 



The quantity es{r,t) is the field emitted per unit area at the source but whose amplitude 
includes implicitly an inverse distance dependence. Propagation is described by the quantity 
g{r,t,i',rs), which is the propagator for a point source at location rg. It includes a phase factor 
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for free-space propagation as well as phase and amplitude factors associated with DISS.Q Here v 
is the center frequency of the passband with bandwidth Az^ that is selected by the receiver and 
mixed to baseband. The notation for e(r, t) and n(r, t) leaves this center frequency implicit; but 
each of these quantities is a time series of a narrowband process and is the baseband equivalent 
of the narrowband radiation field selected by the receiver. Further justification for this model is 
given in Appendix A. 

The field emitted by the source is taken to be amplitude modulated noise (Rickett 1975; 
Cordes 1976b), es(rs,t) = a(rs, t)m(rs, t), which is a physically motivated and empirically 
confirmed model for most astrophysical sources. The amplitude modulated noise model includes 
nonstationary modulations a(rs,t) of stationary noise m(rs,t). The noise correlation function is 
(m(rs]^, ti)m*(rs2, ^2)) = '^(rs2 ~ rsi)A(t2 — ^i), where the asterisk denotes complex conjugation, 
angular brackets denote ensemble average, and A(t) is a continuous delta-function-like quantity 
with A(0) = 1 and width equal to the reciprocal bandwidth of the receiving system. The additive 
noise, n(r, i), also is 'A' correlated in time and is assumed spatially uncorrelated across nonzero 
baselines, b = r2 — ri. 



3.1. Interferometer Visibility Function & Phase Structure Function 

It is well known that the mean visibility function of a scattered point source (in strong 
scattering) is the product of the true source visibility rs(b) and the second moment of the DISS 
modulation, 7g (e.g. Rickett 1990 and references therein): 

(r,(b,ri)) = (e(r,t)e*(r + b,f + ri)) = A(ri)7g(b,0,0,0)r,(b). (3) 

The source visibility is the usual Fourier transform of the brightness distribution, 
Is{ys) = A(rs) = a2(rs), 

r,(b)= j (ir,e+^'='^"'^-^/,(rs), (4) 

where d is the source-observer distance and k = We use a spatial vector, Fg, to define 

the source brightness distribution rather than using an angular variable, as is common practice. 
We assume that the propagation delay between the pair of sites has already been removed, so the 
visibility function maximizes at Tj ^ (Ai^)^"*^. 



^ The form of Eq. |l| is approximate. Factoring the integrand relies on Es being narrowband with bandwidth 
<^ V, that it varies much faster than the DISS propagator, g, and that the bandwidth is much smaller than the 
characteristic scintillation bandwidth, i.e. the characteristic bandwidth on which g changes. This last constraint 
is identical to requiring that differential propagation times be much less than the shortest characteristic time scale 
of the signal. Later, we also consider variations of g with frequency. The signal model presented still applies if we 
consider the total frequency range to comprise many separate intervals in each of which g is piecewise constant. 
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The factor 7g in Eq. |3| is proportional to the second cross moment of the propagator, 
g(r, t, z/, Fg), at two observation positions, times and frequencies separated by b, r, 6iy, respectively, 
and for two point sources separated by drg-. 

rg{h,T,6iy,5r^) = {g{r,t,iy,r^)g*{r + h,t + T,iy + Su,rs + Sr^)) = e-'^-fg{h,T,5u,5r^). (5) 

In this equation, the phase '0 is determined by free-space propagation and drops out in much of 
what appears below, but is responsible for the Fourier relation in Eq. ^. Appendix ^ shows the 
details. The form for 7g with zero frequency lag {di^ = 0) is simply expressed in the Gaussian limit 
using the phase structure function, Z?<^(b, r, (Jrg), 

jg{h,T,6u = 0,6r^) = e-^^^^^'^'^^'^l (6) 

For nonzero frequency lags, a closed-form expression for 7g is not usually available. The DISS 
"gain" G = \g\'^, has unit mean, (G) = 1, and has a normalized autocovariance in the strong 
scattering (Rayleigh) limit, 

7G(b, r, 6iy, 6r^) = {G{r, t, u, r^)G{r + h,t + T,iy + 5iy,r^ + 6r^)) - 1 = |7g(b, r, 6i^, 6r^)\'^. (7) 

For a medium in which scattering occurs with variable strength all along the line of sight, the 
phase structure function is (Lotova &: Chashei 1981; Cordes & Rickett 1998) 

D^{h,T,6r,) oc [''dsC^is)\Ks{sT (8) 
Jo 

bcff(s) = {s/d)h + V^s{s)T + {l-s/d)6rs (9) 

Vefr(s) = (s/(i)Vobs + (l-s/(i)Vp-V„,(s). (10) 

is the coefficient in the wavenumber spectrum for electron density variations and a is the 
exponent of the structure function. For a square-law structure function a = 2, while for a 
Kolmogorov medium in strong, but not superstrong scattering (Cordes &: Lazio 1991), a = 5/3. 
Vp is the pulsar velocity, Vobs is the observer's velocity and Vm is the velocity of the scattering 
material in the ISM. For the case G^{s) oc 6{s — Dg), where Dg is the distance of a scattering 
screen from a pulsar, we retrieve the result applicable for a thin-screen. 

Note that the phase structure function oc for radio propagation through tenuous plasmas. 
For optical and infrared (IR) propagation through the atmosphere, oc A~^. 

Our expressions Eq. 0-|l^ are quite general, being based on Gaussian statistics for the 
wavefield and on the saturated (Rayleigh) regime of scattering. As such, they can be used for 
observations of Galactic and extragalactic radio sources, including pulsars, masers, microquasars, 
active-galactic nuclei, and gamma-ray burst sources. 

Our form for applies to isotropic scattering irregularities. Evidence exists for anisotropics 
in heavily scattered sources (Frail et al. 1994; Wilkinson, Narayan & Spencer 1994; Yusef-Zadeh 
1994; Molnar et al. 1995; Desai & Gwinn 1998; Spangler & Cordes 1998; Trotter, Moran & 
Rodriguez 1998). Though, for simplicity, we consider only the isotropic case in this paper, it is a 
simple matter to extend our results to the anisotropic case, which we will do elsewhere. 
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3.2. Isoplanatic Scales 

DISS is correlated over spatial and temporal scales at the observer's location that are 
determined by contours of constant Dfp. For a thin screen at s = Ds from the source, we 
write = (|befj|/6e)", with h^s given by Eq. |^ and where be is the 1/e scale of the structure 
function. We define the isoplanatic length scale and time scale through D0(biso,O,O) = 1 and 
D^{0, Atd,0) = 1, yielding biso = be{d/Ds) and, if the source's speed dominates Veff (as it does 
for many pulsars), Aid = &e/(l — Ds/d)Vp. We also define the isoplanatic scale at the source's 
location using ^(/,(0, 0, (5rs iso) = 1, resulting in drg^iso = &e(l — Ds/d)~^ = VpAt^. The isoplanatic 
scale drg^iso defines the separation at which two point sources would scintillate with a correlation 
coefficient of e^^. This scale determines whether DISS can resolve a source, as discussed in the 
Introduction. For reference, the Fried scale ro, defined in the optical and IR literature, is related 
to our definitions using D^{b,0,0) = 6.88{b/rof/^ (e.g. Goodman 1985), so ro ~ 3.26e- Also, in 
the radio case rg oc X~^^^ while rg oc A^^/^ for optical/IR propagation. 

In strong scattering, the isoplanatic scale is smaller, in some cases by three orders of magnitude 
or more, than the Fresnel scale. The Fresnel scale, at meter wavelengths, ~ V XD ~ 10^^ cm for 
kiloparsec distances. 

The isoplanatic scale 5rs^iso is smaller for sources that are scattered more heavily. This 
corresponds to more distant sources, sources viewed through regions of excess scattering, or 
sources observed at longer wavelengths. For a continuous source, such as one with a Gaussian 
brightness distribution, its size compared to drg^iso determines the depth of modulation of the 
DISS. To apply this basic idea, however, averaging over time and over the receiver bandwidth 
must also be dealt with carefully because it too affects the depth of modulation. We now consider 
all these effects in what follows. 

4. Resolving Sources with Scintillations 

Several methods can be used to exploit the scintillation phenomenon in order to resolve 
sources. These are (1) measurement of the fractional modulation of the source through analysis of 
the intensity variance; (2) estimation of the intensity PDF, equivalent to analysis of all moments; 
and (3) use of cross-correlation functions for the DISS of separate sources to measure the spatial 
offsets of those sources. 

4.1. Visibility and Intensity Statistics 

The visibility function and the intensity are both second field moments. Estimates of these 
second moments from finite data sets fluctuate by amounts that are formally described by the 
fourth field moment. Encoded in these fluctuations is information about source structure and 
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the intervening medium. We use a normalized fourth moment — a generahzed modulation index 
(squared) — to cast intensity and visibility fluctuations in a similar form. The modulation index 
includes the effects of averaging over time and frequency and is calculated for an arbitrary source 
brightness distribution. 



4.2. Autocorrelation Functions 

To model realistic cases, we take into account averaging over time and the finite bandwidth of 
the narrowband signal. For simplicity, we refer to the finite bandwidth as "frequency averaging.' 
To resolve DISS, the averaging intervals T and B must be smaller than the characteristic 
correlation scales of G, the DISS gain in time and frequency. These are usually called the DISS or 
'scintillation' time scale and bandwidth, denoted At^ and Afj, respectively. 

Let I(r, t) be the intensity at location r calculated for a narrowband signal with bandwidth 
B after averaging over the interval [t — T/2,t + T/2]. Define also the averaged visibility function 
as r(b, i), calculated between two sites separated by baseline b and with time lag r but with zero 
frequency lag. The autocorrelation functions (ACFs) of the intensity and visibility are 

i?j(b,T) = (T(r,t)I(r + b,t + r)) (11) 
%(b,r) ^ (r(r,t)r(r + b,t + r)) (12) 



4.3. Modulation Indices 

To quantify fluctuations of I and F, we calculate the modulation indices, 



,„J(l.,.) . (13) 

where we normalize by the mean intensity in both cases. As shown in Appendix ^, the modulation 
index receives contributions from four terms. For compact sources, the dominant term is in fact 
caused by DISS and is given by 

"^iss(b,r) = {ly^ / (irsi(irs2/s(rsi)/s(rs2)<9iss(b,r,rs2 -rsi,r,S), (15) 
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where, for uniform averagingp] in T and B, 



<5lSs(b,r,5rs,r,B) 



(TB)- 






r' 
















) 




T 




'b 





7G(b,r' + r,5z^, Sr^ 



intensity 



(16) 



ikd ^h-Sr^ 



7G(0,r' + r,(5i/, Sr^ 



visibility. 



with 7g defined in Eq. 0. Visibility fluctuations are independent of baseline b for unresolved 
sources, for which the complex exponential 1, while intensity variations depend more strongly 
on b, even for sources unresolved by the baseline. The baseline-independent property for visibility 
fluctuations is similar to the conclusion found by Goodman & Narayan (1989). The difference 
between intensity and visibility statistics arises from the ordering of the time-averaging and 
cross-correlation operations in the two cases. 

The modulation index, as presented here, includes only DISS fluctuations. There are 
additional contributions to the visibility or intensity variance from intrinsic source fluctuations and 
from additive noise. These are secondary to our discussion here, but are important in any practical 
application where the time-bandwidth product is low and where intrinsic source fluctuations are 
high. For pulsars, pulse-to-pulse amplitude variations are important when only a few pulses are 
included in any averaging. Appendix |^ gives full expressions for all contributions to intensity 
variations. 



4.4. Number of Degrees of Freedom in Fluctuations 



TB averaging and extended structure expressed in the integrals of 7g in Eq. |15|,|16| diminish 
scintillation fluctuations. The modulation index of the averaged intensity or visibility depends on 
time averaging and source extension in similar ways because both increase the number of degrees 
of freedom in the integrated intensity. The number of degrees of freedom is 

A^dof = 2mis| = 2A^iss > 2, (17) 

where Aiss is the number of independent DISS fluctuations ("scintles") that are averaged. For 
observations in the single speckle regime, where scintles are resolved in time and frequency, we 
expect 1 < iViss <, 2. 



^Uniform averaging over frequency corresponds to spectrometer passbands that are perfectly rectangular; real 
passbands, hiy), have shapes similar to Gaussian functions and can be handled by replacing 1 — \5i' / B\ in Eq. ^ 
with the autocorrelation of h{v) and extending the limits to zbcxa. Computations show that this refinement produces 
no significant changes in our discussion. 
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4.5. Examples 

In Figures |l[]3| we show ?Tijgg plotted against integration time T for different values of 
bandwidth B. These cases are for a point source. In the figures, we use T and B normalized by 
the scintillation time scale, Aid, and bandwdith, Af^- Figure I is the case for a thin screen with 
a square-law structure function (i.e. a = 2). Figure |2| is for a thin screen with a = 5/3, the form 
appropriate for a Kolmogorov medium if the scattering is strong but not ultrastrong (Cordes & 
Lazio 1991). Finally, Figure ^ is the Kolmogorov case for an extended medium with uniform 
statistics along the line of sight. 

The differences between the plotted curves for the different media are subtle, but nonetheless 
significant if one were to use the modulation index (for given T, B, say) to try to determine what 
kind of medium was responsible for a given measurement. More importantly, the differences must 
be considered when assessing whether the source is extended. 

To consider source-size effects, we adopt a circular Gaussian brightness distribution with 
source size, Ur, 



/,(r,) = {2naf.) exp (^-^ j . (18) 

For this case, the squared modulation index is 

m?ss(b,r) = (47ra2)-iy'd5r,e-(l^'-^l/2-'-)'QlSs(b,r,5rs,r,i?). (19) 



Figures show mf^g plotted against source size in units of the isoplanatic scale, S^s^iso 
(defined in § |3.2| ). The three figures are for the same square-law, thin-screen Kolmogorov, and 
uniform Kolmogorov media considered in Figures |^^, which are results for point sources. The 
different curves are for different combinations of T/Atd and B/Av^- One feature to note is 
that w-jgg falls off more rapidly with source size <Tr/Srs^iso than it does with T/Atd or with 
B/Ai'd- This is because the source-size dependence results from a two-dimensional integration 



over the difference vector — fsi in Eq- 15 as compared with one-dimensional integrals for TB 
averaging. These figures show again that any inference on source size must account not only for 
TB-averaging, but also for the type of medium underlying the measurements. Furthermore, the 
predicted contributions from TB-averaging rely on accurate measurements of the DISS time scale 
and bandwidth. 



4.6. Interpretation of Modulation Index 

Application of Eq. ^ is as follows. If mf^^ = 1 (within errors) then the source is unresolved 
by the DISS, the baseline b has not resolved the scattering disk, and the scintillations cannot have 
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decorrelated over the averaging intervals T and B. The DISS gain G then has an exponential PDF 
associated with the two degrees of freedom in the scattered wavefield. 

Alternatively, m-igg < 1 can signify (1) resolution of the scattering disk by the baseline (in the 
case of intensity interferometry) ; (2) variation of the DISS over the averaging time or averaging 
bandwidth; or that (3) the source has been resolved by the DISS, i.e. that it is comparable to or 
larger than the isoplanatic scale of the DISS. To discriminate between these possibilities, auxiliary 
information is needed that characterizes the dependence of 7g on its four arguments, b, r, 5u^ and 
6rs- Such information is obtained by making DISS and angular broadening measurements over a 
wide range of frequencies (e.g. Rickett 1990). 

Complications in estimating mjgg arise from the fact that scintillating sources fluctuate, on 
inverse-bandwidth time scales and on a variety of longer time scales, and there is additive noise 
in any real-world receiver system. We consider all such complications in Appendices ^ and |^ and 
also in the next few sections. 



4.7. PDF of the Averaged DISS Gain 

While the modulation index of visibility fluctuations may allow inference of source structure, 
an analysis of the full probability density function (PDF) may be more sensitive. Here we 
investigate the PDF for several cases. 

First consider a scintillating point source with no TB averaging of G. As is well known, the 
intensity PDF is a one-sided exponential in the limit of no additive noise because G is a chi-square 
random variable (RV) with two degrees of freedom, (e.g. Goodman 1985). TB averaging and 
extended sources increase the number of degrees of freedom, and therefore decrease mjgg, as 



discussed in § |4.4| . If TB averaging and source superposition are viewed as combining statistically 

^iss' 



independent RV, G is distributed as 'x^p^ 



where r(x) is the gamma function and U{x) is the unit step function. In detail, however, the 
intensity (or visibility) is the integral of variables that are statistically dependent, so the xIa^iss 
PDF is only an approximation to the true PDF of G. 

The true PDF is calculated by solving a homogeneous Fredholm equation of the second kind 
(Press et al. 1992, pp. 779-785) that results from expanding the propagator g'(r, t, rg, i^) onto 
an orthonormal set of eigenvectors (e.g. Goodman 1985, pp. 250-256) and requiring that the 
expansion coefficients be statistically independent (a Karhunen-Loeve expansion). The general 
case, where time-bandwidth averaging and source extension must be considered, requires solution 
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of the eigenvalue problem (see Appendix 

r-t+T/2 r+B 



rt+T/2 r+B f \f)u\\ C i 

(TB)-'J^_^^^ dt' ddu (l-^j J dr^i [/s(rsi)/.(rs2)]'/' 

7g(0, t" - t', 6u, - rsi)lpnit', Vsi) = Xn1pn{t' , rs2), (21) 



where the A„ are the eigenvalues. The time and frequency integrands are slightly different because 
the wave propagator, g, is integrated over frequency before squaring of the wavefield, while 
time-averaging occurs after squaring. 

For the simple example of a time average of G at discrete times tj = jAt, j = 1, . . . ,N, the 
eigenvectors ipn{tj) and eigenvalues A„ are solutions of 

N 

Y.7giO,ti-tj,0,0)Mti) = >^nMtj), n = l,...,N. (22) 
1=1 

The PDF of G = -/V^^ X^i [^(^i)!^ is the A^-fold convolution of one-sided exponential PDFs, 
each of which has mean Xn/N, because each coefficient in the expansion of g{t) is 
statistically independent and is a complex, Gaussian RV. The PDF can be written in 
the form faiG) = J2n^nGwi~G N/Xn) where the c„ are functions of the eigenvalues, 
Cn = X^^ Hn'^ni^ ~ Xn' / Xn)~^ ■ Frequency averaging behaves similarly. 

Goodman (1985) has shown that exact PDFs calculated in this way are fairly well 
approximated by the X2N1SS with the appropriate number of degrees of freedom given by 



Eq. |17|. The approximate PDF has the same mean and variance as the exact PDF. In the limits 
-^iss 1 and A^iss ^ li the two PDFs become identical. Also, if an observable is calculated as 
the sum of strictly independent DISS fluctuations with equal variances, then the exact PDF is the 

xItViss PDF- 

Figure ^ shows exact PDFs for different instances of time averaging and finite source size. 
The results are for a Kolmogorov wavenumber spectrum (a = 5/3) but do not differ substantially 
for a square-law structure function. Figure ^ is a sequence of PDFs for time-averaging and 
for a point source, while Figure |^ is a sequence of PDFs for finite source sizes but with no 
time-averaging. Also shown in Figure |^c are approximate PDFs based on the PDF for different 
numbers of degrees of freedom given by 2A^iss, where values of A'^iss are chosen to yield PDFs of 
similar variance as in Figure |^,b. 

Comparison of the panels in Figure |^ indicates that time averaging and source extent produce 
similar forms for the PDF of G. This conclusion verifies the notion that, from a statistical point 
of view, TB-averaging and source extension produce like effects in scintillation fluctuations. The 
figure also supports the notion that one may use the approximate PDF to make calculations 
rather than solving the Fredholm equation for every case. This is useful because in cases where 
time and frequency averaging as well as source extent are important, the Fredholm solution may 
not be obtainable. 
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4.8. PDF of Visibility Fluctuations 

The PDF for G derived in the previous section excludes contributions from source fluctuations 
such as those that arise from the amplitude modulated noise model. Here we give a nearly exact 
treatment that takes into account all fluctuations. For simplicity, we assume that observation 
baselines do not resolve either the source or the seeing disk from the scattering. In this case, we 
write the average visibility across baseline bjj between two sites i and j as 

r ^ G{1) + {Mi)5i, +X + C, (23) 

where (/) is the mean source intensity, 5ij is the Kronecker delta, {J\f) is the mean of A/i = 
(the background noise intensity), X is a real Gaussian RV with zero mean, and C is a complex 
Gaussian RV with zero mean. Source fluctuations are described by X which includes noise 
fluctuations associated with m(rs,t) and amplitude fluctuations associated with a{vs,t) in the 
amplitude modulated noise model. C includes additive radiometer noise combined with source 
noise fluctuations, but is uninfluenced by source amplitude fluctuations. Expressions for a\ and 
u'q are given in Appendix C.2. 



The PDF for the visibility magnitude is calculated by successively integrating over the PDFs 
for the different, independent terms in Eq. as done in Appendix ^ The PDF for the scaled 
visibility magnitude, 7 = |r|/(Tc is 

/,(7) = f dOfciG) f dXfx{X)Ue-l^^'^'''^')lo{ii)\ (24) 

J J L J *=((I>+-^/<^)/'^C 

where /q is the modified Bessel function. The integrand factor in square brackets is the 
Rice-Nakagami PDF of a signal phasor added to complex noise (e.g. Thompson, Moran & Swenson 
1991, p. 260). 

To demonstrate the importance of various terms and factors, we first show, in Figure |^ the 
visibility PDF when we vary the bandwidth, taking into account that Vjss varies as we do so. 
There is a tradeoff in discerning the underlying shape of the DISS gain PDF (which encourages 
use of narrow bandwidths) and maximizing S/N, which favors larger bandwidths: as we decrease 
the bandwidth there is less averaging of the scintillations but contributions from noise (the X and 
C terms) increase, thus widening the PDF. 



We compare the X and C terms in Eq. 23 with the G(l) term, which dominates when the 
source is strong. First, we calculate the PDF for |r| with various terms excluded. Figure]^ shows 
the visibility PDF when there are no DISS and no intrinsic fluctuations, just additive noise. 
Figure ^ shows the visibility PDF with the intrinsic fluctuations turned on for a fairly high S/N 
observation, but still with no DISS variations. These curves indicate that source fluctuations can 



contribute significantly to the shape of the PDF. Figure 11 shows the PDF for different source 
intensities (top panel) along with (bottom panel) the difference between the true PDF and the 
PDF where pulsar fluctuations are ignored. The difference vanishes when the source intensity is 
zero and ~ 1% for flnite source intensities. The error is largest near the peak of the PDF. 
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Figure 12 shows the PDF when DISS is included with varying numbers of degrees of freedom, 
2A''iss, but with constant bandwidth. This case would apply to observations of sources with 
different intrinsic sizes or to a point source observed with varying amounts of time averaging. As 
-^iss oo, the PDF tends toward a Gaussian form. The results indicate that the net PDF for |r| 
is extremely sensitive to the number of degrees of freedom in the DISS. Also, it is evident that 
a pure point source can have statistics that mimic those given by an extended source if there is 
sufficient TB averaging. 



Comparison of Figures |12| and |g indicates that some of the changes in PDF shape evident in 
Figure ^ that might be due to, say, effects of source size, are indeed masked by any changes in 
resolution bandwidth in the signal processing. 



4.8.1. When is Pulsar Noise Important? 

In § |C.2| we give expressions for different terms in the visibility variance, including those 
involving pulsar fluctuations and one involving radiometer noise. All fluctuations diminish, of 
course, with increased averaging time. However, the relative sizes of the fluctuating terms are 
independent of the averaging time. In studies of the shape of the visibility or intensity PDF, pulsar 
fluctuations will be comparable to radiometer noise fluctuations when the pulsar signal strength 
satisfles 

where (I) is the source flux density and Wj is its characteristic time scale (either intrinsic or 
imposed by some sampling scheme); B is the bandwidth. For a pulsar, Wj would be the sample 
window in pulse phase and (I) the flux density in that window. 



4.9. Cross-Correlation Functions 



As described in Appendix [B.7| , when the intensity is measured for two sources and cross 



correlated, the total cross modulation index includes only an ISS term because noise and source 



fluctuations do not correlate. Defining the cross-correlation of the average intensities as in Eq. B27 
and defining the cross modulation index as 

2 _ 12 ('^i^) ~ -^1-^2 

"^cross(b,r) = ^ — , (26) 

I1I2 

we find that ?7T,^j.Qgg is identically equal to mjgg of Eq. 

The cross modulation index can be used in two ways to detect source structure or spatial 
offsets between two sources. First, if its maximum is smaller than unity, quenching of DISS by 
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source structure is signified. Also, if the cross-modulation maximizes at a time lag r 7^ for zero 
baseline (b = 0), that signifies a significant offset between sources. We show these cases in the 
Appendix ^ and will apply this method to pulsar observations in a later paper. 



5. Inferring Source Size from Intensity & Visibility Statistics 

Given measurements of the visibility of intensity that have time and frequency resolutions that 
resolve DISS, it is possible to place constraints on source size by comparing visibility fluctuations 
with those expected from the time and frequency resolutions and as a function of source structure. 
Here we outline general Bayesian procedures that first analyze visibility /intensity data and then 
are restricted to a simpler analysis of just the modulation index. 



5.1. Likelihood Analysis of Visibility & Intensity Fluctuations 

Given a set of visibility (or intensity) measurements 

{T^{hi),i = l,N} (27) 

we can normalize them by the rms off-source noise, ac, and use the PDF of 7 = T/ac (c.f. Eq. 24) 
to calculate a likelihood function 

^ = n/7(70- (28) 

i 

The likelihood function depends on numerous parameters that could be estimated by maximizing 
C. These include parameters that describe the source, &s, such as source structure and flux 
density; those that describe wave propagation through the ISM, ©ism, including distances, type 
of medium (e.g. Kolmogorov), and scintillation parameters (Ai/^, At^); and those that characterize 
the receiver and telescope system, 0r, including the telescope gain and system temperature. 
Many of these parameters will be known from auxiliary observations. 

Denoting all parameters collectively as = (©s, ©ism, ©r), we identify the data probability 
P(T>\@) as the likelihood function and we calculate, in standard Bayesian fashion (e.g. Gregory & 
Loredo 1992), the posterior PDF for the parameters as 

^ " Jd@P{®)P{V\®)~ Jd®C' ^ > 

where P(©) is the prior PDF for the parameters and the denominator normalizes the PDF. The 
second equality follows if we assume that the parameters have a flat prior PDF. In the case where, 
a priori, we know many of the parameters, we adopt delta function priors and marginalize them 
by integrating over those parameters. For example, if we wish to derive the PDF of only source 
parameters, we would integrate over (i©isM '^©R to obtain P{@q,\T>). 
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To apply this approach we need to solve the multidimensional Fredholm problem that includes 
source extent and time-bandwidth averaging if we want to use the exact PDF for scintillation gain. 
Alternatively, we could use the approximate PDF based on the PDF (Eq. ^0|) by calculating the 
effective number of degrees of freedom associated with source extent and TB averaging. Another, 
simpler, approach is to analyze only the second moment of the visibility /intensity fluctuations, as 
we now consider. 



5.2. Inferring Source Size from Modulation Indices 

The squared modulation index m-jgg can be calculated through direct estimation of moments 
or by fitting a PDF shape to a histogram of visibilities (or average intensities). Such estimates of 
rftjgg are typically made from data that span a large number of scintles, A'iss.total S> 1. This 
number is approximately 

"~»(i + c|^)(i + c|^), (30) 

where the total observing time and bandwidth may be written as Ttot = NtT and i?tot = NyB, 
using T and B as the basic resolutions in time and frequency defined previously. The characteristic 
time and frequency scales for DISS are At^ and Afj, respectively. The factor (" k, 0.2 — 0.3 
takes into account that scintles are not packed tightly in the v — t plane. In Cordes (1986), I 
conservatively used C = 0.1 whereas a more accurate calculation yields the values presented here. 

The fractional estimation error on m-igg due to the finite number of scintles is 
(7^2 K, 2mjggA^jgg^r^Qrp^L- Invoking the Central Limit Theorem for A'iss,totaL) we expect 
TTijgg to have PDF, A^(mjgg, (Tm2). 

We write the likelihood function in terms of the PDF for wijgg estimated from data and using 
the model for ^rijgg (as given by Eq. [Tsj) : 

C = [2710^2 j exp ["iiss(0s) - "iissj J > (31) 

where ©s is a vector of parameters that represents the source structure. Bayes' theorem can then 
be applied according to Eq. |2^ to derive the posterior PDF for 0g. For the Gaussian brightness 
distribution of Eq 18, the posterior PDF is simply a one-dimensional PDF fari'^r) for the sole 
source parameter ar. 

Uncertainties in the application of this inference scheme include the systematic errors 
associated with not knowing the true form of the structure function for the medium and also the 
statistical errors in the measured modulation index, w,jgg and in the DISS parameters Ai/d and 
Aid- We address these uncertainties in Paper II. 
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6. Summary and Conclusions 

We have derived a general methodology for analyzing diffractive interstellar scintillation 
fluctuations that is applicable to single aperture and interferometric observations. In this paper, 
we considered only the strong scattering regime where the scattered wavefield has Gaussian 
statistics. The method explicitly takes into account timc-bandwidth averaging that is often used 
in the statistical analysis of such observations. Such averaging modifies the statistics in a way that 
is identical to the effects of extended source structure. We show that time-bandwidth averaging 
and extended source structure both increase the number of degrees of freedom in the scintillations 
from the minimum value of two that describes the fully modulated, Gaussian wavefield of the 
scintillations. 

Our methodology can be applied to any radio source in the strong scattering regime, including 
compact active galactic nuclei and gamma-ray burst afterglows. In another paper, we will address 
sources of these types and we will also consider scintillations in the weak and transition scattering 
regimes. 

In Paper II we apply our results to the recent VLBI observations of the Vela pulsar by 
Gwinn et al. (1997) and find that the scintillation statistics may be accounted for fully by 
time-bandwidth averaging. Any contribution from extended source structure is less than an upper 
limit of about 400 km at the 95% confidence interval. This upper limit on the transverse extent 
is substantially larger than the size expected from conventional models that place radio emission 
well within the light cylinder of the pulsar and close to the surface of the neutron star. 

I thank Z. Arzoumanian, S. Chatterjee, C. R. Gwinn, H. Lambert, M. McLaughlin, and B. 
J. Rickett for useful discussions and H. Lambert and B. J. Rickett for making available their 
numerically-derived autocovariance functions for Kolmogorov media. This research was supported 
by NSF grant 9819931 to Cornell University and by NAIC, which is managed by Cornell University 
under a cooperative agreement with the NSF. 
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APPENDICES 

A. Scintillating Amplitude Modulated Noise Model 

Here we derive a general statistical model that incorporates incoherent summing in the source 
and wave propagation throughout the interstellar medium. 

The narrowband (scalar) electric field incident on an aperture and selected by a feed antenna 
and by a bandpass receiver may be written in the form 

EA{t) = Re {e{t) exp(-ia;ot)} , (Al) 

where ujq is the center frequency and the complex, baseband wavefield is e (e.g. Thomas 1969). 
The baseband field is often explicitly extracted through quadrature mixing schemes in heterodyned 
radio receivers (e.g. Thompson, Moran &; Swenson 1991, p. 150). 

Early work on pulsars modeled £{t) as amplitude modulated noise (AMN) with additive 
background and receiver noise: 

e{t) = a{t)m{t)+n{t), (A2) 

where m and n are complex, Gaussian wavefields that describe intrinsic source noise and additive 
noise, respectively. The factor a(t) is a real modulation function that describes source variations 
on time scales much longer than the reciprocal center frequency or reciprocal bandwidth; 

otherwise the statistics of a{t) are arbitrary. The additive noise has a modulation index 

m7v = [(|m|4)/(|m|2)2-l]i/2 = i. 

The AMN model was first presented by Rickett (1975), who attributed complex, Gaussian 
statistics to m{t). For this case the modulation index is also unity, tum = 1- Cordes (1976b) 
considered Poissonian shot-noise statistics for m{t) based on physical models for pulsar emission; 
for Poissonian noise, mM > 1- 

Empirical tests on pulsars (Cordes 1976a; Hankins Sz Boriakoff 1978; Gordes &; Hankins 1979; 
Bartel & Hankins 1982) show consistency of m{t) with Gaussian statistics on time scales as short 
as ~ 1 /xs. Tests in the time domain resort to investigation of the relative amplitudes of various 
terms in the autocorrelation function of the intensity (Rickett 1975; Cordes 1976a; Bartel & 
Hankins 1982). Tests in the frequency domain, with the same conclusion, use the autocorrelation 
function of the spectrum (Cordes & Hankins 1979). The model predicts that there is frequency 
structure in the spectrum of a single pulse with characteristic bandwidth equal to the reciprocal of 
the time duration of a{t). This frequency structure averages out as multiple pulses are summed. 

Tests on OH and H2O masers (Evans et al. 1972; Moran 1981) show that maser emission 
also conforms to the AMN picture. It is expected that any radio source can be described by AMN 
because large numbers of particles contribute to the observed signals and most, if not all, natural 
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sources involve incoherent superposition of radiation from incoherent or coherent emissions from 
individual radiators. Thus AMN should apply to gamma-ray burst sources and the most compact 
AGNs that show intra-day variability. 

A.l. Amplitude Modulated Noise for Extended Sources 

We model extended sources as follows. First, the baseband field produced by a point source 
at location (rs,z = 0) is 

£s(rs, t) = a(rs, t)m{rs, t), (A3) 

where a(rs,t) is the (real) amplitude modulation and m(rs,t) is complex Gaussian noise (Rickett 
1975). Here and everywhere, vectors are two dimensional and perpendicular to the line of sight. 
The quantity is the field emitted per unit area at the source, uninfluenced by propagation 
(either through free space or a turbulent medium), except that we include the dependence on 
distance from the source on its mean amplitude, for simplicity. For a steady source, a{rs,t), is 
constant in time. For pulsars, it describes the periodic envelope of pulses that modulates the 
underlying noise process, m{rs,t). The corresponding mean intensity is, using A = a^, 

Is{r„t) = (|£.(r„i)|') = {A{v„t)). (A4) 
In most of this paper we assume stationary statistics, so Is{rs,t) — /s(rs)- 
The total measured field is the integral over source components 

£(r, *) = / *s £s(rs, t)g{r, t, v, rg), (A5) 

where we include a multiplicative propagation factor, (/(r, t, zv, Fg), defined in the next section. 
The noise, with stationary statistics, has a correlation function 

r2m(rsi,rs2,ii,*2) = {m{rsi,ti)m*{rs2,t2)) = 5(rsi - rs2)A(t2 - h), (A6) 

where the asterisk denotes conjugation and A(r) is an Hermitian function having unit amplitude, 
A(0) = 1, and width approximately equal to the inverse of the receiver bandwidth. Angular 
brackets denote an ensemble average, except where noted. The noise fourth moment is the 
standard dual sum of products for a complex Gaussian process, 

r4m(rsi, rs2, rs3, rs4, ti,t2, h, U) = (m(rsi, ti)m*(rs2, t2)w-(rs3, *3)"^*(rs4, ^4)) 

= r2m(rsl, rs2, il, t2)r2m(rs3) rs4, ti) 

+ r2m(rsi,rs4,ti,i4)r2„(rs2,rs3,i2,i3)- (A7) 
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A. 2. Propagation Through a Thin Diffracting Screen 

Consider the following geometry: a point source at (rs,0), a thin screen at {r',Ds) and an 
observer at {r,d = Dg + D). The screen changes only the phase of incident waves. Under the 
narrowband approximation (Az^ ^ u) (so that all phase factors may be considered constant over 
the band) the propagated baseband field for a point source at rg is (e.g. Goodman 1985) 

e(r, t, rs) = (iXD)-' J dr' e'*^""' ''-'"'^^^h'^'^'^e,{T,, t - c-'Vn), (AS) 

where D = {Dg^^ + D^^)^^ and the integral is normalized so that a screen with zero phase yields 
simply a delayed version of the emitted field, es{TCs,t — c~^d). Under the paraxial approximation 
(transverse scales much smaller than line-of-sight distances), 

Vis = V12 + V23 



2D 

V23 » D+ 



r — r'P 



2D 

We assume further that variations in propagation times, c~^Pi3 and c~^I?23 (as a function of 
relevant source locations rg and screen exit points r') are negligible compared to the characteristic 
variation time scales for (j){r',t) and e{rs,t). c~^Az/(Pi3 — d) <C 1. Typically, Eg varies on time 
scales of order the reciprocal bandwidth (e.g. 100 fis or less) while g varies, due to the changing 
geometry, on time scales of seconds to hours or more, for the situations we wish to consider. 
Therefore, for Ss to be factored out of the integral, we require Av{T)i3 — d) <C 1. This simply 
means that any time smearing from differential arrival times must be less than the time resolution 
of the signal. Though we assume in the remainder that the inequality is satisfied, we point out 
that there are many instances where it is not, corresponding to the well known 'pulse broadening' 
effect (e.g. Rickett 1990). As a rule of thumb, when pulse broadening is important, scintillations 
are difficult to resolve in time and frequency. And when scintillations are important, the pulse 
broadening can be too small to be important, as we assume here. 

We can now write the propagated field for a single point source as 

e(r, t, rg) = es{rs,t - d/c)g{r, t- D/c, u, r^), (A9) 



where g is the propagator. 



g{r,t,i/,rs) = (iXD) ^ J dr' exp ^^i 



^(l?.-^|r'-r,|2 + D-i|r-rf)+0(r',t) 



(AlO) 



and we use k = 27:c~^v. The normalization of g yields it to be simply a unit modulus phase factor 
when the screen phase is zero and also it yields (l^l^) = 1 when the screen phase has Gaussian 
statistics. In the following we will ignore the delays d/c and D/c in Eq. in our notation. 
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The time dependence of g{r,t,i',rs) arises from motions of source, observer and medium, which 
influence all terms in the exponent in the integrand. Absent any random phase screen (i.e. ^ = 0), 
g is simply a complex phase factor that describes free-space propagation, 

g{r, t, V, rs) = e^'=l'---'-|'/2d. (All) 



A. 3. Propagator Second Moment for a Thin Screen 

The propagator's second moment across a baseline b, at two times separated by r for two 
point sources at 2 and at two frequencies separated by is 

Tg{h,T,6v,5Ts) = {g{Y,t,v,Ys)g*{r + h,t + T,v + 5v,rs + 5ys)) 

= e-^^7g(b,r,(5z/,5rs) (A12) 

^ = A:d-H(rs-r)-AR+^|AR|2} (A13) 

AR = [(Jrg-b + (Vp-Vobs)r]. (A14) 

There is no term oc Au in ijj because it is negligible according to the narrowband assumption made 
earlier, i.e. that c~^Ai/{Vis — d) ^ 1. 

The real, second moment 7g for zero frequency lag is 

7g(b,r,5i/ = 0,5rs) = e-5^^(b'^'^"--) (A15) 
where the phase structure function and its arguments are given by 

D^{h,T,dr,) ^ ([<^(r,0-</.(r + beff,t)f) = (A16) 

beff = {Ds/d)h + V,sr+{D/d)6r^ (A17) 
Veff = {Ds/d)Vohs + il-Ds/d)Vp-Vm (A18) 

The scaling exponent for the structure function, D^, is a, which takes on a value a = 5/3 for a 
Kolmogorov spectrum in some instances. The e"^ scale of \jg\^ is be, Vp is the pulsar velocity, 
Vobs is the observer's velocity and Vm is the velocity of the scattering material in the ISM. The 
normalized autocovariance function for G = \g\'^, which we call the intensity "gain," is (in the 
strong scattering, or Rayleigh, regime) 

7G(b,r,(5i/,(5rs) = (G(r,t, z/,rsi) G(r + b, t + r, i/ + rsa) - 1 = |7g(b,r,(5i/,(5rs)p, (A19) 

and the mean intensity is 

7(r, t) = J dvs Is{rs, t) = j drs (^(rs, t)). (A20) 
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Note that {G) = 1. We emphasize that the form for 7g in the second equahties of Eq. A12 and 



Eq. A19| rehes on the assumption that the scattered wavefield is sufficiently randomized that 



Gaussian statistics apply. 

For nonzero frequency lags, 7g generally must be obtained through appropriate numerical 
integration (e.g. Lambert & Rickett 1999; Lee & Jokipii 1975). For the special case of a square-law 
structure function, a closed form expression is available for a thin screen (Chashei & Shishov 1976; 
Cordes et al. 1986; Gupta et al. 1994). 



A. 4. Extension to An Arbitrarily Thick Medium 

The results presented so far were derived explicitly for a thin screen. For a medium in which 
scattering occurs with variable strength all along the line of sight, the results may be extrapolated 
quite simply and generally. Assuming that the net measured field is still Gaussian, it is a matter 
of simple geometry to work out the form of the equivalent phase structure function. Following 



Lotova &: Chashei (1981) and Cordes & Rickett (1998) we use the same results as in Eq. A12 - A26 
but make the replacements, 

D^{h,r,5v,) = {Xreffa dsCl{s)\hesisT (A21) 

Jo 

beff(s) = {s/d)h + V,s{s)T + {l-s/d)6r, (A22) 

Veff(s) = (s/d)Vobs + (l-s/d)Vp-V^(s), (A23) 

where 

•' a2"r(l + a/2) ^ ^ 

and is the coefficient in the wavenumber spectrum for electron density variations (Cordes & 
Lazio 1991; Armstrong, Rickett & Spangler 1995). For the case C^(s) oc 5{s — Dg), we retrieve the 
thin-screen results of the previous section. 



A. 5. Visibility Function 

The ensemble mean visibility function is the product of the true source visibility and the 
propagator's second moment (ignoring a phase factor) 

{Te(h,t,Ti)) = (e(r,t)e*(r + b,t + ri) = A(ri)7g(b,0,0,0)r,(b,t), (A25) 

where we have designated the interferometer lag as r^. This must match any geometrical time 
delays to within the reciprocal of the receiver bandwidth. The source visibility is the usual Fourier 
transform of the brightness distribution, 

r,(b,t) = I drse+*'=^"'^-V.(rs,i). (A26) 
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The time lag Tj in Eq. A25 is zero to within a very small time (of order the reciprocal bandwidth; 
see Thompson, Moran & Swenson 1991). By contrast, the time lag in Eq. A12 extends over very 
long times, seconds to hours, that characterize DISS fluctuations. 
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B. Intensity & Visibility Fluctuations 

Here we consider variations in the time-averaged intensity and visibility for the model 
of Appendix We derive the modulation fractions of these quantities taking into account 
scintillations and intrinsic and additive noise. Our expressions will include any averaging of the 
scintillation modulation over frequency as well as time. 

We define the general fourth moment for the narrowband field e [sans additive noise) 

Riei^MMMM) = (e(r,tiK(r,f2)e(r + b,t3)e*(r + b,f4)>, (Bl) 



where e(r, i) is given by Eq. A5. Expanding out, i?4e involves fourth-order moments of a(rs,t) and 
m(rs, t) from the AMN model of Appendix ^ and of the propagator g{T,t,v,Ts). We use the fact 
that a,m, and g are statistically independent. Also, m and g are complex gaussian processes so 



their fourth moments are dual sums of products of their second moments, as in Eq. A7 



The fourth moment of of a(rs,t) ends up as the second moment of its square, A = a?, which 
we assume has the form 

(^(rsi,ti)A(r,2,t2)) = (A(rsi,ti))(A(rs2,t2))[l+m^/0A(rs2 -rsi,t2-ti)]. (B2) 

This form assumes stationary statistics for A with a correlation function pA and modulation index 
niA- As we show in paper II, assuming stationary statistics is not restrictive, even for pulsars 
which are highly nonstationary across pulse phase but appear to have stationary statistics when a 
fixed pulse phase is considered. For other radio sources with intrinsic variations much longer than 
those of pulsars, we may consider {A[rs,t)) to be constant in time, with mA = 0, at least over a 
typical observation time of minutes. 

Our treatment also includes any variations of the scintillation propagator, 5, across the 
bandwidth of the narrowband signal, e(r, i). By partitioning e into subbands in which the 
propagator is piecewise constant and between which the emitted signal £s is statistically 
independent, we easily can incorporate finite bandpass effects while being consistent with our 
earlier assumption about the narrowband signal.^ 



B.l. Autocorrelation Functions of Time Average Quantities 

The time-averaged intensity, 

/■t+T/2 

I(r,t) =r"i / dt'l{r,t'), (B3) 

Jt-T/2 



^The net effect is that intensities from the subbands add while the intensity autocorrelation function discussed 
later involves an average over frequency lag. 
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has autocorrelation function 



i+T/2 



i?j(b,r) = (I(r,t)I(r + b,i + r)) = T'' dtadtbRi{h,tb - ta + t) (B4) 

J Jt-T/2 

= r-i|_^^dr'(^l-^^Rj(b,r' + T). (B5) 

The integrand is given by 

Rl{h,T) = Rie{h,h,h,ti + T,ti + t). (B6) 

Similarly, we consider the visibility, 

r,{h,t,Ti) =e{r,t)e*{r + h,t + Ti) (B7) 

and its time average, 

_ , rt+T/2 _ 

T{h,t,Ti) = T-^ dt'T{h,t',Ti). (B8) 

Jt-T/2 

The autocorrelation, analogous to Eq. B5, is 

i?r(b, r; r,) = T'^ J^^ dr' (^1 - ^) R^{h, t' + r; n), (B9) 

and involves the integrand 

i?p(b, r;r) = i?4e(b, ti, ti + r, ti + r + r, + r). (BIO) 

Note that we distinguish here between the lag associated with the definition of the visibility, Tj, 
and the lag r with which we consider the autocorrelation of the visibility. 



B.2. Modulation Indices 

We are most interested in the normalized variances of the time-average intensity and visibility. 
These are defined in terms of the autocorrelation functions as 

mj(b,r) = ^^^2 (Bll) 



2 _ i?r(b,r)-|(r,)p 
m^(b,r) = , (B12) 



where we normalize by the mean intensity in both cases. 

The total modulation index squared for the intensity or visibility is the sum of three main 
terms: 

m^(b,r) = mi^ss(b,r) +mpsR(b,r),+m^OlSE(b,r), (B13) 
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where ?7iigg measures the contribution from scintillations only, rripo,^ measures the contribution 
from source amplitude fluctuations combined with scintillations, and Jtinoise measures source 
noise fluctuations. Later (§ |B.6| ) we will also consider the effects of additive radiometer fluctuations. 
For pulsars, ^pgj^ includes pulse shape variations and noise fluctuations. For sources that are 
steady over an observation span of minutes to hours (or more), Jripg^ = 0. The 'noise' term, 
"^NOISE' depends on source structure and corresponds to the output of an intensity interferometer 
that is proportional to the square of the visibility function (e.g. Hanbury-Brown 1974, pp. 48-49). 

We can write these terms as 

mhsi^,r) = (I)-2||dxdyIr,(x)Ir,(y)Qiss(b,r,y-x,r,i?) (B14) 

mi.sR(b,r) = (I)-2y'y'dx(iyIr,(x)Ir,(y)QpsR(b,r,y-x,T,i3) (B15) 

(b,r) = (I)-'y'y'dxdyIr,(x)Ir,(y)QNOiSE(b,r,y-x,r,B) (B16) 



m 



NOISE 



The 'Q' functions are defined in terms of integrals over time lag, like that in Eq. and over 
similar frequency-lag integrals, that we denote as 

{X{y))y,Y = dy - M j x{y). (B17) 



B.3. Intensity Fluctuations 

For intensity fluctuations, the Q functions are 

(3lSs(b,r,5rs,r,B) = 

(7G(b, r' + r, 6iy, 5rs))^,_y. ^^^^ , 
QpsR{h,r,6r„T,B) = 

{(PA('5rs,r' + r))^,^y + (pA((^rs, + r)7G(b, r' + r, (5rs))^,^^. 



(B18) 
(B19) 



5u,B 



+Ra{t, T)e 
QNOiSE(b,r, 6rs,T,B 



ikd ^b-5rs 



'PA 



{5v,, 0) [7G(b, 0, 0, 0) + (7g(0, 0, 5v, 5v,))^^^j,] } , 



-ikd ^b (5rs 



= i?A(r,r)e 

We have made use of the lag-integrated noise correlation 

Ra{t,T)^T-^ f^dr' (l-^)|A( 



7G(b, 0, 0, 0) + (7g(0, 0, 6i^, dr,))^^ ^ . (B20) 



-'m2 



(B21) 



Recall that A(r) is a function with unit maximum amplitude [A(0) = 1] and width of order the 
reciprocal bandwidth, B^^. For integration times T ^ B^^, we have 



RA{r,T) 



Wa 
T 



UiT-T)U{T + T), 



(B22) 
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where U{x) is the unit step function and W/\ = J dr |A(r)p ~ i? is the characteristic time scale 
of the noise fluctuations. 

Note that for T ^ Wa, the terms involving i?A in QpsR and Qnoise are much smaller than 
Qiss- Also, for radio sources other than pulsars, Qpsr = because = 0. 



B.4. Intensity Interferometry 

We can relate our results to those of Hanbury-Brown and Twiss (e.g. Hanbury Brown 1974), 
who used the intensity autocorrelation function to determine the magnitude of the source visibility 
function for optical stars. The term of interest in their work corresponds to our Qnoise, in 
particular the second term which involves (7g(0, 0, 61^, 6rs)) ^ ■ The first term with 7G(b, 0, 0, 0) 
vanishes for baselines much larger than the Fried scale. Also, the apertures used by Hanbury 
Brown and Twiss were larger than the Fried scale, causing aperture averaging that we have 
not treated but which is analogous to time-bandwidth averaging. For ground-based optical 
observations of stars, Hanbury-Brown and Twiss used a bandwidth such that scintillations were 
constant over the band and the stars they observed were much smaller than the isoplanatic scale. 
In this case, 7g(0, 0, (5rs) 1 and the effect was maximized. Note also, however, that the 
amplitude of the effect scales with R^{t,T) Wa/T (BT)^^, which is small for significant 
time-bandwidth averaging. 



B.5. Visibility Fluctuations 

For the visibility we have, using AV = Vp — Vobsj 
Qiss(b, r, 6r,,T, B) = | A(r,)|' (70(0, r' + r, 5v, 5^^y-^kd-^^-[^r.+^^(r'+-) 

'3pSR(b,r,(5rs,r,S) = 

m\ {1 A(r,)|2e-*^'^"'b-^'-= {pA{5v,y + r)7G(b, 0, 0, 0))^, ^ 
+ {pa{6v.. r' + t)7g(0, t' + r, s^^y-ikd~^^-[^w{r'+r)] 

+RAir,T)p{Sr,,0) [l + {jG{h,n,6u,6r,))s,^B]} ■ 
QNOlSE(b,r,5rs,r,S) = Ra{t,T) 1 + (7G(b, r^, ^rs))^^^ 



t',T; 5i^,B 



,(^23) 



(B24) 



t'T; Su,B 



(B25) 



B.6. Effects of Additive Noise 

Results given so far for intensity and visibility fluctuations have considered only the signal 
emitted by the source. Including the additive radiometer noise n as in Eq. we obtain an 
additional contribution to the total modulation index, that we denote nif.^^: For the time-average 
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intensity and visibility, the contribution is 

{2{l)~'^{Ni)5ijRA(T, T) (l + i(I)"^(Ni)) intensity fluctuations 
(B26) 
(I) " ^ (N, ) i^A (r , T) ( 1 + + (I) " ^ (Nj ) ) visibility fluctuations 

We have labelled the baseline with ij indices to represent the ith and jth sites. The Kronecker 
delta indicates that the contribution for the intensity holds only for single-site measurements for 
which i and j are equal. 



B.7. Cross Correlations of Time- Average Intensities 

In some circumstances, we are interested in the cross correlation function of the intensity 
between two sources that may or may not be scintillating together. Pulsars, for example, have 
different pulse components that may come from spatially different emission regions and it is 
possible to record or compute intensities for each component separately and cross-correlate them. 
We define the cross correlation as 

q^,2(b,r) = (Ti(r,t)l2(r + b,t + r)). (B27) 

The utility of the cross correlation is that it is affected by both the separations of the two sources 
and the size of each source. For example, the time lag at which the CCF maximizes is determined 
by the separation of the sources and by the effective velocity, Veg (Eq. 10|jAlg| ) . 



The cross correlation simplifies greatly if we assume that the amplitude-modulated noise 
in each source is statistically independent from the other. For pulsars, this is a reasonable 
assumption in many cases, though in others where there are drifting subpulse fluctuations that 
appear successively in different pulse components, this is an approximation. Letting Ii^2 be the 
(ensemble) mean intensity of each source, we find that 

Ci^^(b,r) ^ h{t)h{t + T) (B28) 



For two point sources, one at Vgi, another at rs2, that have stationary statistics, the 
normalized crosscovariance is 

CjJh,T)-hl2 
hh 



l-i,, (b, r) ^ ' = (7G(b, r' + r, Sv, r,2 - r^i)),, ^ . (B29) 



To illustrate the utility of the crosscovariance consider the case where 7g is constant over the 
averaging intervals T and B. The time lag that maximizes 7j^^ is the solution of 

^L»4b,r' + r,rs2-rsi) =0. (B30) 
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For a thin screen this becomes (c.f. Eq. |6| 

D 



d_ 
dr 



which has the solution 



Veff • [[Ds/d)h + {D/d)5v, 

Veff' 



(B31) 



(B32) 



The abihty to estimate Tmax with precision depends on its value relative to the characteristic DISS 
time, which is the width of 70 as a function of r. Defining Atj using 7g(0, Atj, 0, 0) = f 
find that Atj = ^d/Veg, where Id is the characteristic diffraction scale, yielding 



we 



Aid 



Veff • [{Ds/d)h + {D/d)5v, 

4Veff 



(B33) 



For b = 0, we expect to identify Jrg 7^ only if {D/d)5rs is a sizable fraction of id and also if 
the effective velocity is not orthogonal to Srg. The definition of 'sizable' depends on the number 
of independent ISS fluctuations used in any estimate of the cross correlation function, which is 
A^iss given by Eq. |l^. The error on Tmax ~ ^t^Niss'^^'^ , so a three-sigma measurement requires 
Srs\\ ^ miss-'^\did/D), where = ■ Veff. 



Smirnova, Shishov &; Malofeev (1996) give a similar expression for Tn 
medium with a square-law structure function. Our result is more general. 



that is based on a 
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C. Probability Densities for Strong Scattering 

Here we derive probability density functions for tiie DISS gain in strong scattering. We use 
an exact treatment based on Karliunen-Loeve expansions that take into account source extent and 
time-bandwidth averaging. Then we derive the PDF for the visibility and intensity that takes into 
account all features of the amplitude modulated noise model of Appendix ^ 



C.l. Exact Solution for the PDF of the Scintillation Gain 

The time average intensity may be written as 

T(r,t) =T-i / dt' f dr,h{r,)\g{r,t',u,r,)\^ (CI) 
Jt±T/2 J 

for a source with arbitrary brightness distribution Ig modulated by DISS. The DISS modulation g 
has been bandwidth averaged in accord with the considerations of Appendix ^ 

g{v,t,u,r,)=B-^ ( du' g{r,t,u' ,v,). (C2) 

Jv±B/2 

1/2 _ 

Following the approach described by Goodman (1985; pp. 250-252), we expand Is (rs)^(r, i, i/, rg) 
onto a set of orthonormal basis vectors ipnit^Ys) with coefficients hn- The orthonormality condition 
is 

/ /dr,Vn(t,rs)C'(i,rs) = <5„n' (C3) 

Jt±T/2 J 

and the 6„ are given by 

bn=T-^ I dt [ dr,ll/\r,)g{r,t,u,r,)Mt,rs). (C4) 
Jt±T/2 J 

By requiring that = (|6nP)<Jnn' (i-e- that the 5„ are statistically independent), the following 

eigenvalue problem results: 



(TB)-' J^_^^^ dt' j^^ d6u (l-^j J dr^i [/.(rsi)/s(rs2)]'/' 

7g(0, t" - t', 5u, rs2 - rsi)i^n{t', r^i) = A„Vn(i', rs2), (C5) 



where An = (l^nP) are the eigenvalues. The time and frequency averaging are handled differently 
because the wave propagator, g, is integrated over frequency before squaring of the wavefield, 
while time-averaging occurs after squaring. 

The expansion implies that 

I(r,t)=^|6„P (C6) 
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and that 

(/(r,t))=^A„. (C7) 

n 

The expansion coefficients are gaussian distributed because the integral Eq. C4 is a sum of 
gaussian variables. Therefore, each term in Eq. |6nP) is exponentially distributed and the 
intensity PDF is the convolution of each of these exponentials. 

The convolution can be calculated through Fourier transforms and inverted using the residue 
theorem yielding, for nondegenerate eigenvalues, 

N 

fi{I) = Y.Cne-'/^-U{I), (C8) 

n=l 

where U{I) is the unit step function and the coefficients are given by 

TV 

Cn = Xn-^ n (1 - V/A„)~' . (C9) 

When only one eigenvalue is important, as it is for a point source with negligible time-bandwidth 
averaging, the PDF for I contains only a single term with mean A = (I). 



C.2. Visibility PDF for the AMN Model 

Here we present an alternative derivation of the visibility PDF that takes into account all 
source, propagation and additive-noise fluctuations. 

The visibility function is the product of the narrowband fields from two sites (i and j)0, 

T{t)=ei{t)e*{t), (CIO) 

where = gijam + Uij for a point source that produces an identical field at the two sites. The 
propagator and the additive noise are both different at the two sites, in general. 

The instantaneous value of the visibility is 

r = giQjAM + nin* + a{gimn* + g*m*ni). (Cll) 

We have simplified the notation, using i,j to label spatial location rather than using location and 
baseline vectors as we have used in previous sections. The ffist term is the scintillated pulsar 
signal, the second is due to additive noise at the two sites, while the third term represents cross 



*We take the product at the same time in order to keep notation simple. In practice, a delay must be introduced to 
account for the different optical path lengths to the two sites. Our notation assumes this has already been removed. 
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products. If there were no scintillations, the visibility would be a noisy phasor (from the pulsar) 
combined with complex noise. Scintillations modify the source phasor to make it complex, in 
general. However, the pulsar noise (from A and M) are in phase with respect to the scintillations. 

The (ensemble-average) mean visibility is 

{T) = {gig*){A) + {Ni)6ij. (C12) 



In practice, a time average is used to approximate the ensemble average, with attendant 
errors. We use the following notation for the time-averaged visibility: 

, rt+T/2 

T = {r{t))BT^T-^ Tit), (C13) 

Jt-T/2 

where the subscript 'BT' on the angular brackets denotes time averaging of a bandlimited process 
with bandwidth B. 

By expanding A and M into mean values and zero-mean fluctuations, e.g. A = (A) + 5A and 
M = (M) + 5M, we can write the time-average visibility as 

r = {g,g*) ^^{A){M) + {Ni)5,j +X + C, (C14) 

where 5ij is the Kronecker delta. The first term is due to the source, the second term is the mean 
system noise for a single site observation {i = j), and the last two terms are fluctuations, 

X = {g^g*[5A + 6M{A) + 6A5M])BT 

C = {a{gimn* + g*m*ni))BT + {riin*) - {Ni)5ij. 

We separate X and C because in useful limiting cases, discussed in the next subsections, they 
become, respectively, real and complex processes. For i = j, C also becomes real. Moreover, X is 
in phase with the phasor term, {gigj)^^{A){M), while C is randomly phased. 

If there are no intrinsic fluctuations, X vanishes and C then depends only on additive noise. 
However, the AMN model demands that there be source fluctuations even if there are no amplitude 
modulations. We let (M) = = 1 without any loss of generality. 



Though we have assumed a point source to arrive at Eq. pl4 , the equation also applies 



to extended sources that are spatially incoherent. Spatial incoherence yields summation of 
contributions from different source elements that also imply Gaussian statistics. 

All terms in X and C are uncor related, so variances of individual terms sum to yield the 
total variance. We assume time-bandwidth averaging such that BT ^ 1; thus X and C become 
Gaussian random variables (GRVs) by the Central Limit Theorem. However, we assume BT is 
small enough so that the DISS factor, {gigj) ^j,, is not a GRV. For now, holding gi{t) and gj{t) as 
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fixed realizations of the DISS fluctuation (i.e. not averaging over an ensemble for these quantities), 
we find that X has PDF 0, A^(0, crx^)) while C has PDF Nc{Q^ac'^), where the variances are 

ax^ = {\X\'') = {k)^{G,Gj))j,^{BT)-^[BWAm\ + ml,{l + m\)\ (C15) 
^c' = (|C|2) = (2i3r)-i[(A)((N,)(G,)^^ + (N,)(G,)5^) + (N,)(N,-)], (C16) 

and we have used Gi = etc. The forms of these variances are consistent with expressions 
given by Rickett (1975). In general, we can write the variances of the real and imaginary parts 
of X and G as ax^^ = ^^^(l ± Pdj) and cr^^ . = |cr^(l ± 5ij). When the DISS is perfectly 
correlated between the two sites, the correlation coefficient pc^ , (which is equal to 7G(b, 0,0, 0), 



c.f. Eq. A19 ) is unity and X is real. As the DISS decorrelates between the sites, the — > crxr- 
Only when the two sites are identical (e.g. for a single aperture measurement of intensity) is G 
real. For all interferometers, G is complex with equal variances of the real and imaginary parts. 

Deriving the variances involves assumptions about the correlation times for the various signal 
terms and how they influence mean squares of the time averages. We have assumed that n and m, 
the noise processes, have correlation times much smaller than that of the amplitude modulation, a, 
which in turn has a much smaller correlation time than the integration time used. This hierarchy is 
consistent with the fact that the noise correlation times are the reciprocal of the bandwidth used. 
For pulsars, data are often obtained by using only a small range of pulse phase but averaging over 
many pulse periods. Pulsar pulses are broad band but decorrelate on times about equal to the spin 
period. Our expression for ax uses the correlation time Wa for the amplitude modulation. This 
is effectively the width of the correlation function pA defined in Eq. |B^. As applied in Paper II, 
we would take Wa to be the width. At = PA(pp, of the pulse window used in the analysis, where 
P is the pulse phase and A(j)p is the window width in pulse phase units (cycles). Then T = NpAt. 
We have also included niM, the modulation index of pulsar noise fluctuations. For amplitude 
modulated noise with Gaussian statistics, niM = 1- By retaining it, we can see what changes in 
the statistics if we artifically turn off the noise fluctuations, ac depends on (M) but not on mj\,/. 

Note that the contribution of pulsar noise to visibility fluctuations relative to the contribution 
from additive noise is independent of the averaging time. Consider, for example, the ratio 
ax/o'+, where a+ is the value of ac when there is no source. We have ax/cr+ oc SNRq, where 
SNRq = (A) / ^ (Nj) (Nj) is the ratio of source strength to system temperature, when both are 
in the same units (e.g. Janskys). Similarly, ac/cr^ — 1 oc VSNRq. This result is at odds with 
Gwinn et al. (2000), who state that pulsar "self noise" can be ignored because it diminishes with 
averaging time. It does diminish but it cannot be ignored unless the signal to noise ratio is small. 

The PDF for F given gi{t) and gj{t) may be calculated by appropriate integration over the 
Gaussian PDFs for X and G. We now consider some specific cases. 



^A^(0, ax^) denotes a Gaussian PDF of a real variable with zero mean and variance ax^ while Nc{0, ac^) denotes 
a complex Gaussian quantity having real and imaginary parts with equal variances, ac^. 
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C.3. DISS Perfectly Correlated Between Sites and Constant over BT 

For perfectly correlated DISS (between sites i and j), {gig*) {Gi) ^j, = {Gj) ^j, = {G) ^j, 
and (GiGj)) {G"^) ^^p^ thus X becomes real. If, moreover, the DISS modulation is constant 
over the averaging time T, then {G) ^j, G = constant and (G'^)^y G^ . The visibility becomes 

V = G{k) + {^i)5,j+X + C. (C17) 

We assume a point source (and strong, saturated scintillations) so that, with G constant over T, 
the scintillation PDF is a one-sided exponential. The PDFs of individual elements of F are 

/g(G) = e-'^UiG) 

fx{X) = N{0,ax^) 

fc{C) = iVe(0,ac') 

ax^ = G^{A)\BT)-^[BWAml + ml,{l + m\) 

ac' = (2i3r)-M((N.)(N,)) + +G(A)((N,) + (N,))]. (C18) 



C.3.1. PDF of the Complex Visibility 

Let T = Tr + iTi. The PDF of F is 
/r(r) = J dGfciG) J dXfx{X)fc{Tr - G{A) - X,Ti) 



1 



(F, - G(A) - Xf + f'' 



= J dGfciG) J dX(27rax)"^/^e-^'/24(27rac)-^e ^^^c ^ (dig) 
Performing the integral over X, we obtain 

MF) = (2wc)-^e-^?/-^ / dGfa{G)ial + alr^^e-l^^' " ^^^^^'/^^^^ + ^^).(C20) 

Note that for no signal ((A) ^0), we get 

/r(r) = (2vr4)-ie-|riV24, (C21) 

a circular Gaussian PDF. Our expression in Eq. |C20| disagrees with Eq. 11 of Gwinn et al. (2000), 
which assigns equal variances to the real and imaginary parts of F. The variances are not equal, 
in general. Also, there is an extra factor of 2-k in their equation. 



0.3. 2. PDF of the Visibility Magnitude 
The PDF of |F| can be calculated as 

/|r|(|r|) = |r| r"#/r(|r| cos ,/.,|r| sine/.). (C22) 
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Here we take a slightly different approach. It is convenient to scale the magnitude of the 
visibility by the rms of the complex term, ac- Using 7 = |r|/crc and i = {A)/ac, the conditional 
PDF for constant G and X = is 

^(7|i, G) = 7e-^(^'+G^^^)lo(7i), (C23) 

where Iq is the modified Bessel function. This result is the well known Rice-Nakagami PDF for 
a real phasor of length i combined with a complex Gaussian phasor (e.g. Thomson, Moran &: 
Swenson 1991, Eq. 9.37). Integrating over the PDF for X, we have, for fixed G, 

f^{j\G) = J dXfx{X)f,{j\i + -^,G). (C24) 

Then, integrating over the PDF for G, we have the PDF for 7 that takes into account all 
fluctuations, including DISS, 



^(7) = J dGfGiG)f^{^\G). (C25) 

For S/N -^00, ac^ 0, 

F = G{A) + {Ni)6ij + X, (C26) 
and the PDF of |r| for fixed G becomes N{G{A) + (Ni)5ij,ax'^), with the PDF for 7 given by 



Eq. C25 



C.4. Perfectly Correlated DISS but G 7^ Constant over BT 

Specializing to the case of a weak source for which ax ^ crc, we have (using G = (G)^j.) 

F « G{A) + {Ni)5ij + G 

^""^ ' GT{N,ss) ^ 

where F(x) is the gamma function and U{x) is the unit step function. G is distributed 
approximately as X2N ' ^ chi-square random variable with 2A'igs degrees of freedom, where 



^iss 

A'lss = TTT-iss s-nd where mjgg is given by Eq. 17. The true PDF of G is obtained by solving the 



appropriate Fredholm equation for the eigenvalues that determine the PDF, as described in the 
main text. 
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Table 1: Symbols and Acronyms Used 



Symbols 


Definition 


/. . .) 


KlncipTn nip iWTf^vi^ ctp 


(■ ■ ')bt 


T'lmp fivprj^pfp nvPT" t.imp T of nrnpp*?*? wit, In VifiTiHwiHt.Vi R 

-L line Ct V d&jC w V d UllilC -L yji. CKi yJ^- \J\j\j\j\j W 1 lilX UCtliVJ. VV iXJi uH -I—/ 


ACF 


A iThnporrplatinn fnri ptinn 


CCF 


OrossroTTpl ation fnn rtion 


DISS 


Diffractive Interstellar Scintillation 


DM 


Disnprsion Mpasiirp 

J — ' A.VJ L/ \^ J. kJX VJX-L XV X V_><^kJ I^X x_> 


RISS 


Rpfrartivp Tntprstpllar Srintillatinn 

X li^^XX CXiV>' UX V ^ XXX U^X O UVjXXCtiX V^'XXXuXXXCXiUXVyXX 


n(f) 


A nTnli'l'nriP TyiAnnlj^tiAn 1 tpj^ 1 i 


A = a^ 


TnfiPTi^^ifv monn Ifitiinn 

XlllJt;iXOX \j\ LLWJyJi. Lli-CtUlV^lJ. 


B, Au 


R p Ti n WT^ H f n 


b 


Rpcplirip 
C(jO X X X X 


h 


1 /p cfalp of nViJ^^p ^t.rn ft.n TP fnnptinn 


baft 


effective baseline 




Intensity CCF 




OopfficiPTit in plprtroTi-HpTiRitv wfivPTimnliPr snPctrnTn 

Vw/V_/V/XXXV^X\^XX(J XXX V^X"^*^ UX WXX VJ.V/XXkjX U y VV Cv V '-/XXLA.XXXI^/ V/X OL/^-'^UX LLXXX 




T^j^rt ii-smirpp Hmt^iTiPP 

1 It li I uXX OVJ IJ.X VjO 1_XX0 LOjXX\^C 


D = d — Dg 


Q p-ppp-ri _ T?. Q T*-J- n n 1 cf a rt fp 


^ s 


Smirpp-Rcrppn Histfmcp 




Phase structure function 


J a 


Numerical factor in structure function 


fr(G) 


PDF of scintillation modulation 


f (^) 


PDF OT nnrmi^li'/prl viRinilitiv 

J. J y A- \JL XXVJX XXXCIjXXZJt^VJ. VXOXkyXXXLjy 


o 

i) 




G = lop 


Scintillation modnlation of intpnsitv or scintillation "srain " 

vyxxx uxxxcti uxv^xx xxxv^ i_cxcti uxv^xx vyx xxx u v^xxux u y vyx lj v^xxxuxxxca uxv^xx c^ca>xxx ■ 


Ea 


IVarrowliaTirl plpptrip fiplH 

X I CtiX X \J VV k^CtiXXU, v_^xv_^ ux x^^- xx^^xu. 


£ 


Complex baseband electric field 




Complex electric field at source 


I 


Intensity 


SI 


Intensity flucutation 


I 


Time-average intensity 


jC 


Likelihood function 




Characteristic spatial scale in diffraction pattern 


m{t) 


Complex gaussian noise with unit mean square 


M = |mp 


Squared magnitude of m 


ruA 


Modulation index of A 
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Table 1: Symbols and Acronyms Used (continued) 





T) PTi m ^" 1 r»n 


^ISS 


oClIlTjliicnjlOIl IllOQUlctTjlOIl IIIQGX — rillb / inecLIl j 


ruM 


iViUQUictllUll lIlClcX Ui IVl 


-'*dof 


IN LllliUcI Ui Llc^Iccb Ui iiccLlUiii iii bL/iiibiiidliUiio 




i>l liiilUci Ui iiiUt;pt;iiUt;iil bCiiiUiiidLiUii icdbUico diVcidi^cU. 




A ri (H 1 "1" 1 '\TPi rr^TTi Tr\ 1 p^'v n r^i go 
xA.vJ.U.1 1)1 VC l^Ulli picA. \^ClU.Ol!5lClii iiUlOC 


Tie 


i?ict/ ClcL/UiUil U-cilolLj 


R- 


iT.V>'X' Ui dVci d^cU. iiilciiolLj 


R— 

-£ Lp 


A (ITT r\T Q^rOTti crofi AriGi rii li+Ar 

ii-V^ J- kJL dVCl d^CLl V loiUilit V 


^4e 


± ULli Lli liiUliiClit Ui l^UiiilJltJA. iitJiU 


i s 


T^'xxro ri im on Gi r^n £1 1 '\r(^r*'i' rw sit g/^iitpo 

J- WU U-liliClloiUlidi VCl^LUl dt OULIH^C 


' 




s 


ijUCdliUii diUiig, iiilc Ui Dl^iil), o — U dl oUUiCt; 


T 


Tn'l'OO'T'tti'ir^n rw tnroyfl crm cr I'lmo 
j-iitcgi dijiuii ui dvci dgiiig, Liiiic 


Ti 


± iiiic id^ libcU iii OdiCLiidUiiig, iiitciiciUiiictci Vi&iUiiiUico 


V a 
* cff 






Observer's velocity 


T/ 


r UlSdi VCiULlLy 


T/ 


Vciuciiy UI scdXT/ei iiig niccnuiii 


^SS 


veiucity UI ioo amrdCtiun pattern 


ix 


il/Ji-IJUiiciiU iii JJiidoc oUiUdjUic iliiiCUiUii 


J- I 


1 n't" on Gi'i'TT' Q 1 1 rr^TTQ t*i ci ti i^o tit n /^"J"! <~in 
illUullol L y db lA UU L/U Vdl IdiiL^tJ i Lliilj UiUii 


i £ 


V IbiUlilLy iUliCLlUli 
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oeeuiiu iiiuiiiciiu Ui jjiujjd^diui y 


7G 


dUtUCUVdridllCc UI bClIltllidtlUIl gdlll KJf 




l\r riTm qI 1 voH A liTT r\T nniGO 




l~)iTTTar'f iriTi rvT "^PiTi1"illflf ioti n^iTinwintn 




Diffraction or scintillation time scale 


</> 


Phase perturbation from refractive index perturbations 




Eigenfunction in Karhunen-Loeve problem 


PA 


Autocorrelation function of A{t) 


biso 


Isoplanatic scale of diffraction pattern at observer's location 




Isoplanatic scale at source's location 


(7r 


Length scale in Gaussian brightness distribution 


^iso 


Isoplanatic angle 


ds 


Source angular size 
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Fig. 1. — The DISS modulation index plotted against averaging time in units of the characteristic 
diffraction time, T/At^ for values of the frequency resolution in units of the characteristic diffraction 
bandwidth, B/^v^ = 0,0.5,1,2,5,10. The right hand scale gives the effective number of ISS 
fluctuations that are averaged, A^igg = l/mfgg. This case applies to a thin-screen having a square- 
law structure function. 
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Fig. 2. — The DISS modulation index plotted against averaging time in units of the characteristic 
diffraction time, T/At^ for values of the frequency resolution B in units of the characteristic 
diffraction bandwidth, B/Au^ = 0,0.5,1,2,4,8, as labelled. This case is for a thin screen with a 
Kolmogorov wavenumber spectrum. The right-hand scale indicates the number of ISS fluctuations 
('scintles') averaged, Njss- 
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Fig. 3. — The DISS modulation index plotted against averaging time in units of the characteristic 
diffraction time, T/At^ for values of the frequency resolution B in units of the characteristic 
diffraction bandwidth, B/Au^ = 0,0.5,1,2,4,8, as labelled. This case is for a uniform medium 
with a Kolmogorov wavenumber spectrum. The right-hand scale indicates the number of ISS 
fluctuations ('scintles') averaged, A^igg. 
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T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




J I I I I I I I I I I I I I I I I I I I I I I I I I I L 



1 2 3 4 5 

Scaled Source Size cr / 6r . 



Fig. 4. — The DISS modulation index plotted against source size in units of the isoplanatic scale, 
Srs^iso- The different curves are for different values of averaging time T and bandwidth B in units 
of the characteristic diffraction time scale and bandwidth. In order of light to heavy lines, the 
curves are {T/AU, B/Au^) = (0, 0), (0.4, 0.4), (0.8, 0.8), (1.2, 1.2), (1.6, 1.6), (2.0, 2.0). This case is 
for a thin screen with a square-law structure function. The right-hand scale indicates the number 
of ISS fluctuations ('scintles') averaged, A^iss- 
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Fig. 5. — The DISS modulation index plotted against source size in units of the isoplanatic scale, 
Srs^iso- The different curves are for different values of averaging time T and bandwidth B in units 
of the characteristic diffraction time scale and bandwidth. In order of light to heavy lines, the 
curves are (T / Atd, B / Ai^^) = (0, 0), (0.4, 0.4), (0.8, 0.8), (1.2, 1.2), (1.6, 1.6), (2.0, 2.0). This case is 
for a thin screen with a Kolmogorov (a = 5/3) structure function. The right-hand scale indicates 
the number of ISS fluctuations ('scintles') averaged, A'iss- 
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Fig. 6. — The DISS modulation index plotted against source size in units of the isoplanatic scale, 
Srs^iso- The different curves are for different values of averaging time T and bandwidth B in units 
of the characteristic diffraction time scale and bandwidth. In order of light to heavy lines, the 
curves are (T/Atd, B/Au^) = (0, 0), (0.4, 0.4), (0.8, 0.8), (1.2, 1.2), (1.6, 1.6), (2.0, 2.0). This case is 
for a uniformly extended medium with a Kolmogorov (a = 5/3) structure function. The right-hand 
scale indicates the number of ISS fluctuations ('scintles') averaged, A^iss- 
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0.5 1 1.5 2 

Scintillation gain G 

Fig. 7. — Probability density functions for the DISS 'gain' G for different amounts of time averaging 
and source extent, (a) PDFs vs. averaging time. The different curves going from thinnest to thickest 
hncs are for T/At^^ = 0.01, 0.5, 1, 2, 3, 5, 10 and 20. The PDFs were determined by solving a 
homogeneous Fredholm equation, as discussed in the text, (b) PDFs for Gaussian brightness 
distributions with different sizes relative to the isoplanatic size. The different curves going from 
thinnest to thickest lines are for ar/Svs^iso = 0.01, 0.05, 0.1, 0.2, 0.3, 0.5, 1.0 and 2. These curves 
were determined by solving a two-dimensional Fredholm equation, (c) Approximate PDFs given 
by a X2Afiss '^^^'^ number of degrees of freedom 2Niss, where A^iss = 1, 1-05, 1.125, 1.25, 2, 3, 5 
and 12. 
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Fig. 8. — Probability density function for the visibility magnitude for different bandwidths for 
cases where (A)/(N) = 0.14. The different curves are for different bandwidths, as labelled, which 
correspond to different numbers of ISS fluctuations, also as labelled. This figure corresponds to 
cases where ISS variations and the complex term X are both included in Eq. |C2^^C25| . 
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Fig. 9. — Same as Figure ^ except that ISS variations and the pulsar-noise term have been turned 
off by making the PDF's for G and X delta functions, 5{G - 1) and 5{X), in Eq. p2i - |C25 . 
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Fig. 10. — Same as Figures except that here ISS variations are turned off but the pulsar noise 
variations X are turned on. 
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2 4 6 8 10 

|r| 

Fig. 11. — Top: Visibility PDF for different source strengths. The ratio of signal strength to 
mean noise strength is shown. The PDF's were calculated for Niss = 1-22 and a bandwidth of 25 
kHz. Bottom: Difference between the true PDF and the PDF while ignoring pulsar fluctuations. 
A positive value means that the true PDF exceeds the latter PDF. 
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Fig. 12. — The PDF of the visibihty magnitude when DISS is included with varying numbers of 
degrees of freedom, 2A^iss, for A^iss = 1,1.125,1.25,1.5,2,5,10,20,50 and 100, with some vahies 
labelled. Results are shown for {A) / ^J NiNj = 0.14. As A^iss the PDF tends toward a 

Gaussian function. 



